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ABSTRACT 

Recent observations indicate that dark matter haloes have flat central 
density proflles. Cosmological simulations with non-baryonic dark matter 
predict however self-similar haloes with central density cusps. This contradiction 
has lead to the conclusion that dark matter must be baryonic. Here it is shown 
that the dark matter haloes of dwarf spiral galaxies represent a one-parameter 
family with self-similar density proflles. The observed global halo parameters 
are coupled with each other through simple scaling relations which can be 
explained by the standard cold dark matter model if one assumes that all the 
haloes formed from density fluctuations with the same primordial amplitude. 
We flnd that the flnite central halo densities correlate with the other global 
parameters. This result rules out scenarios where the flat halo cores formed 
subsequently through violent dynamical processes in the baryonic component. 
These cores instead provide important information on the origin and nature of 
dark matter in dwarf galaxies. 

Subject headings: Dark matter — galaxies: spiral — galaxies: structure 
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1. 



Introduction 



The nature of the dark matter (DM) is one of the most important and still unsolved 
problems in astronomy. Current cosmological models assume that there exists a nonbaryonic, 
cold dark matter (CDM) component which consists of non-relativistic particles that are 
assumed to be coUisionless and to interact with the baryons only through gravitational 
forces. CDM simulations indeed reproduce quite well the observed clustering on galactic 
scales (Davis et al. 1985). They also predict the existence of extended dark matter haloes 
around galaxies (Frenk et al. 1985) with constant outer rotation curves, in agreement with 
the observations (Casertano & van Gorkom 1991). 

In order for the galactic gas and stars to achieve a constant circular velocity in the 
outer regions, the density profile of the DM halo must fall approximately proportional 
to in the relevant radius range. Such a profile resembles the density structure of an 
isothermal, self-gravitating system of particles which is characterized by a constant velocity 
dispersion. DM haloes are therefore often approximated by a so called modified, isothermal 
density profile (Begeman et al. 1991) 



where Vc is the core radius and po is the central dark matter density. 

Early cosmological N-body calculations did not have enough resolution to resolve the 
inner structure of DM haloes. Recent CDM-simulations with much higher spatial resolution 
(Dubinski & Carlberg 1991, Navarro et al. 1995) however lead to DM haloes with an inner 
density distribution that is inconsistent with equation 1. Whereas the isothermal density 
profile becomes almost constant at radii smaller than Vc and has a finite central density, 
the numerical models indicate a density distribution which diverges as ~ in the 
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inner parts, leading to an infinite density in the center. Unfortunately such central density 
cusps are hard to verify in normal spiral galaxies as their inner parts are gravitationally 
dominated by baryons. The inferred DM properties then depend strongly on the assumed 
mass-to-light ratios of the galactic disk and of the spheroidal central stellar component. 
The situation is different in dwarf spiral systems which recently have been observed with 
high resolution. It has been shown that some of these systems are completely dominated 
by dark matter (Carignan & Freeman 1988). They therefore represent ideal candidates for 
an investigation of the inner structure of dark matter halos. Flores & Primack (1994) and 
Moore (1994) recently showed that the DDO galaxies have rotation curves which rule out 
singular halo profiles with infinite central densities and are in good agreement with the 
modified isothermal density law (eq. 1). The origin of these shallow DM cores represents 
an interesting and challenging cosmological puzzle which is not yet solved. 



2. The Universal Density Profiles of Dark Matter Haloes 



Scale free cosmological scenarios with non-baryonic, dissipationless DM predict that 
all DM haloes should have similar structure (Navarro et al. 1995). This prediction is tested 
in Figure 1 which shows the DM mass distribution as derived from the observed neutral 
hydrogen rotation curves of four well studied dwarf spiral galaxies. These systems are 
completely DM dominated, even in their innermost regions. Their rotation curves therefore 
trace directly the DM mass distribution and the uncertain gravitational contribution of the 
baryonic component can be neglected. All four profiles indeed follow the same universal 
mass relation, implying a universal halo structure. Note that the modified, isothermal profile 
(dashed curve) only fits the inner regions; it rises somewhat too fast at large radii. The 
dotted and dot-dashed curves show the mass distribution as predicted from cosmological 
CDM simulations on cluster scales (Navarro et al. 1985). These calculations predict DM 
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profiles which are given by the fitting function p[r) = 300 r2oo X r~^[ r + 0.25 r2oo where 
p is the mean density of the universe at the epoch when the DM halo formed. r2oo denotes 
the radius of the DM sphere of mean overdensity of 200. These profiles predict too much 
mass at small radii as a result of the central density cusp (Flores & Primack 1994, Moore 



The observed universal mass profiles can be fitted nicely over the whole observed radius 
range by the following, purely phenomenological density distribution: 



which is shown in Figure 1 by the solid line, po and Tq are free parameters which represent 
the central density and a scale radius, respectively. This revised density law resembles an 
isothermal profile in the inner regions (r ^ Tq) and predicts a finite central density po. Vq 
is similar to the core radius Vc of the isothermal fit formula (equation 1). For large radii, 
the mass distribution of the isothermal profile would diverge proportional to r. Fquation 2 
leads to mass profiles which diverge logarithmically with increasing r, in agreement with 
the predictions of cosmological CDM calculations (Navarro et al. 1995). 



Supposing that the DM haloes of dwarf spirals are indeed all similar with a density 
structure as described by equation 2, we next can investigate the question whether there 
exists a relation between po and Tq. Instead of studying /Jq, which cannot be observed 
directly, we choose as free parameter the rotational velocity Vq of the DM rotation curve 
at Tq. To and Vq are determined by fitting the theoretically predicted rotation curve (using 
equation 2) to the observed DM rotation curves. Accurate neutral hydrogen rotation 



1994). 





3. The Dark Matter Scaling Relations 
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curves are now available for about a dozen low-mass disk galaxies. From this sample we 
select those candidates which are clearly dark matter dominated at Tq. This restriction 
considerably reduces the error in determining Vq due to uncertainties in the mass-to-light 
ratios of the gaseous and stellar components. In addition, the inner DM profiles are not 
strongly affected by the gravitational force of the baryons. Only those galaxies can be used 
for which the rotation curve is known at Tq. These constraints reduce the sample to seven 
galaxies. The present sample consists of the four previously studied galaxies plus three 
systems (Puche & Carignan I99I, Broeils 1992): NGC55, NGC300, and NGCI560. For 
these galaxies the gravitational force of the visible components is not completely negligible 
inside Tq. The contribution of the visible matter to the gravitational potential therefore had 
to be subtracted, leading to a residual DM rotation curve. Here we adopt the best three 
component fits to the rotation curves as suggested by Puche & Carignan (I99I) for NGC55 



Figure 2 shows Vq of all galaxies as a function of Tq. We find a strong, linear correlation 
between Tq and v^'^. Assuming spherical symmetry, one can derive the following scaling 
relations for Uq, the total dark matter mass Mq inside Tq, the central density po and Tq: 



The equations 3 demonstrate, that the density profiles of the dark matter haloes in 
dwarf galaxies depend only on one free parameter, e.g. Tq. Given Tq, the central density po 
(equation 3) and by this the complete DM density profile (equation 2) is specified. One 
should note, that Moore (1994) also noted a trend of increasing size of the galaxy with the 



and NGC300 and by Broeils (1992) for NGCI560. 




(3) 
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peak rotational velocity v^ax- Unfortunately, v^ax is not well determined in most systems 
as the rotation curves tend to rise past the observed region. Moore's conclusion is however 
in qualitative agreement with the present results. 

4. Discussion of the Scaling Relations 

From the scaling laws (equation 3) one can derive the following relation: Mq ~ Vq'^. 
This relation is similar to the TuUy-Fisher- Relation (TFR) for massive spiral galaxies 
(TuUy & Fisher 1977) which predicts that the galaxies' luminosity L correlates with their 
maximum rotational velocity v^ax according to Z ~ v^^^. Note however that the TFR 
combines the mass of the luminous galactic components with the maximum of the rotation 
curve whereas here, the velocity of the DM rotation curve at Tq is correlated with the DM 
mass inside Tq. There exists no obvious reason why the relations should be similar. In 
addition, it has been shown that the TFR is not valid for DM-dominated DDO galaxies 
(Carignan & Freeman 1988) as in these systems the baryonic component is not massive 
enough in order to determine v^ax- 

Like the TFR, the equations 3 can be used as a tool for determining the distances to 
DM-dominated dwarf galaxies, provided their HFrotation curves are known up to Tq. In 
this case, Tq in arcsec can be specified by fitting the universal velocity profile as given by 
equation 2 to the observed rotation curve. Vq can then be determined directly from the 
data. Fquation 3 then gives the corresponding scale radius Tq in units of kiloparsecs. A 
comparison with the observed radius in arcsec leads to the distance of the galaxy. 

There might exist a cosmological explanation for the scaling relations (equation 3). 
Consider a spherical region with DM density p and overdensity 6 = p /]) — 1 in an otherwise 
uniform Finstein-de Sitter universe with mean density ]). The mass of this DM halo is 
Mfi = 4:TrpR^ /3 where R is its physical radius. The density fiuctuation first expands with 
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the universe, reaches a maximum radius at a certain redshift and then collapses 
again. For a coUisionless system of DM particles, this DM halo will achieve a virial 
equilibrium state with a final radius = 0.5i?m- Bardeen et al. (1986) and White (1994) 
demonstrate in detail how one can calculate R^ and function of for certain 

cosmological scenarios. As an example, the solid line in Fig. 3a shows the theoretically 
predicted mass-radius relation for virialized, dissipationless CDM haloes which formed from 
primordial Icr density fiuctuations in the standard CDM-universe. 

Let us now assume that all the observed DM haloes result from Icr fiuctuations. As all 
haloes are self-similar there should exist a universal, linear relation between the observed 
scale parameters (Afo,ro) and the virial parameters (Mh^Rh). In order to determine this 
relation, the dotted curve in Fig. 3a shows the observed mass profile M(r) for the DM halo 
of the dwarf galaxy DD0154. Its scale parameters (Mq, Tq) are shown by the filled circle. If 
DD0154 formed from a Icr fiuctuation, its virial parameters should lie somewhere on the 
solid line. The virial parameters are therefore determined by the intersection (open circle) 
of the dashed and solid curve, leading to the following relations 

Rh = 3.4 X ro (4) 
Mh = 5.8 X Mo 

Supposing that the DM haloes of all the observed dwarf spirals are similar and result 
from Icr CDM fiuctuations, their cosmological virial parameters and R^ can be derived 
from the observed scale parameters Tq and Mq, using the equations 4. Figure 3b shows 
versus R^ for all the dwarf spirals of our sample. The data points indeed follow very nicely 
the theoretical mass-radius relation. 
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5. Conclusions 



As demonstrated in the previous section, cosmological models with non-baryonic, 
dissipationless cold dark matter can in principle explain the scaling relations of DM haloes 
around dwarf galaxies. The dotted and dashed lines in Figure 3b show however clearly 
that a much larger spread of the virial parameters is expected as regions with different 
primordial overdensities could in principle form low-mass DM haloes. If the proposed 
cosmological interpretation is valid we have to draw the conclusion that only fluctuations 
with a certain, fixed amplitude managed to form dwarf spiral galaxies with dominating DM 
haloes. This interesting question will be addressed in greater details in a subsequent paper. 

Cosmological models predict self-similar haloes, in agreement with the present 
observations. Figure I shows however that numerical CDM simulations lead to DM density 
profiles with shapes that disagree with the observations. A possible answer to this problem 
might be dynamical processes within the baryonic, visible component (Flores & Primack 
1994, Moore 1994) which affect the gravitational potential in the inner regions of DM 
haloes, by this leading to a secular change of the DM density profile. Fjection of gas might 
for example decrease the gravitational force in the inner region, triggering an expansion 
of the inner DM halo and leading to a finite central density with an almost constant 
density core. Although such a scenario seems at first very attractive, the present results 
indicate that significant fine tuning would be required between the early cosmological, 
coUisionless collapse phase which leads to DM haloes in the first place and the secular 
energetic processes in the baryonic matter which occur later. It is unlikely that DM haloes 
have self-similar density profiles, if their inner structure is subsequently changed by these 
completely different, secular processes, whereas the outer regions result from the coUisionless 
relaxation phase of the DM halo. We rather would expect to find a second independent 
parameter which describes the structure of the cores of DM haloes and is determined by 
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the efficiency with which the baryonic processes affect the inner DM halo. 

The observed shallow central density profiles of DM haloes are more likely a direct 
result of their coUisionless formation history. If DM, for example, has a finite maximum 
phase space density fmax (as in the case of warm or hot dark matter) conservation of fmax 
during the coUisionless halo formation phase limits the halo's core phase space density /o, 
leading to a central region with /o ~ fmax- If the central velocity dispersion ctq of the DM 
halo is finite, this constraint implies a finite maximum central density po fmax X Cq- This 
scenario would however require that the central phase space densities of DM haloes are 
similar which can be ruled out from the scaling relations. According to the equations 3, /o 
scales as /o ~ Po/'^o ~ f^o^^^ ^^'^ therefore decreases steeply with increasing scale length. A 
finite maximum phase space density therefore is unlikely to be the correct explanation. 

The nature of DM and the origin of DM haloes with finite central densities remains ill 
understood. The discussed discrepancies between CDM calculations and the observations 
indicate that some important physical features, which are related to the nature and origin 
of dark matter are still missing in cosmological models. 

I would like to thank Simon White for interesting and critical discussions and the 
referee for clarifying comments. 
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FIGURE CAPTIONS 



FIG. 1: Dark matter mass profiles, derived from the observed rotation curves are 
shown for the following dwarf spiral galaxies: DD0154 (open triangle, Carignan & Beaulieu 
1989), DDO105 (open square, Schramm 1992), NGC3109 (open circle, Broeils 1990) and 
DDO170 (starred. Fake et al. 1990). The scale radius Tq is defined by equation 2. The 
errorbar at the innermost triangle represents the observational uncertainty in determining 
the rotational velocity and the mass at these radii. For larger radii, the errorbars are as 
large as the size of the symbols. The isothermal fit is shown as dashed curve, assuming a 
core radius Vc = Tq. The solid line shows the revised profile (equation 2). The dotted and 
dot-dashed curves show DM profiles as predicted from CDM calculations. The best fit to 
the data (dotted line) is achieved if one assumes r2oo = 17.5 X Tq which for an = 1, /i = 0.5 
CDM universe corresponds to a formation redshift of z = 0.6. The dot-dashed curve 
shows the mass profile as predicted from the CDM calculations if r2oo = 5 X Tc with a 
corresponding formation redshift of z = 1.5. 

FIG. 2: The scaling relation between the rotational velocity Vq of the DM haloes at 
To is shown as function of the scale radius Tq (equation 2). The open circles show the four 
best known DDO galaxies which have been used in Figure 1. The filled circles show three 
additional galaxies for which the contribution of the baryonic component to the rotation 
curve at Tq had to be subtracted. The error bars show the observational uncertainty in 
determining Vq as quoted in the literature. The dashed line shows a fit (equation 3) through 
the data points. An extension of the data to larger radii fits nicely the halo parameters of 
Sc-Im galaxies (Kormendy & Sanders 1992). For comparison, the dotted curve shows a 
linear fit ~ Tq. 
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FIG. 3: The solid line in Figure 3a shows the predicted mass-radius relation Mh{Rh) 
for virialized, cold dark matter haloes which formed from primordial Icr density fluctuations 
in the standard (0 = 1; h = 0.5)-universe with a biasing factor of 6 = 0.6 (White et al. 
1993). The dotted line shows the halo mass profile for DD0154 as derived from its rotation 
curve and extended to larger radii, using equation 2. The filled circle indicates the location 
of the observed halo scale parameters (Afo,ro), the open circle shows the predicted virial 
parameters(Af/i, for DD0154. 

Figure 3b compares the virial parameters versus of the observed DM haloes 
(open circles) with the mass-radius relation as predicted for a standard CDM universe. 
The dotted, solid and dashed lines represent primordial 0.5(7, Icr and 2(7 perturbations, 
respectively. 



- 16 - 



800 



1 I I I I — \ — I — \ — I — \ — I — \ — I — \ — I — \ — I — \ — r 



I — I 



600 



§ 400 



-X 



O: 



^ o 
> 



200 







-r I \ \ I \ I \ I \ I \ I \ I \ I \ L 

2 4 6 8 10 

[kpc] 



Fig. 2.- 



- 17 - 




log R [kpc] log Rj, [kpc] 



Fig. 3.- 



